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ABSTRACT 

A  formula  for  the  initial  buckling  loads  for  clamped, 
shallow  sphercial  shells  under  uniform  external  pressure 
is  obtained  by  combining  the  solutions  of  two  linearized 
versions  of  the  original  nonlinear  problem.   One  of  these 
versions  is  a  lineeir  eigenvalue  problem,  while  the  other 
is  the  bending  problem  for  a  shallow  cap  in  the  linear 
theory  of  elasticity.   The  formula,  which  is  obtained 
in  a  simple  manner,  yields  buckling  loads  that  are  in 
better  agreement  with  experiments  than  previous  approxi- 
mate solutions  to  the  nonlinear  problem. 
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SYMBOLS 

E  =  Young's  Modulus 

2h  =  Shell  thickness 

k  =  2/3  (1  -  V^) 

p  =  External  pressure,  positive  Inwards 

p  =  Initial  buckling  pressure 

^cr  ° 

P  =  Loading  parameter,  see  equation  [i|.] 

P  =  Value  of  loading  parameter  corresponding  to  p 

R  =  Radius  of  middle  surface  of  shell 

V  =  Vertical  displacement 

w  =  Radial  component  of  displacement 

y  =  Reduced  slope  of  middle  surface  =  a  -  9 

a  =  Slope  of  tangent  plane  to  deformed  surface 

with  the  horizontal  plane 

a  =  a/A. 

n  =  Y  +  P  © 

V  =  Stress   function,    see  equation   [2] 

9  =  Polar  angle 

©  =  ©/TV 

-A-  =  Semi-angle   of  shell  opening 

V  =  Poisson's   ratio 

p  =  Geometric  parameter,  see  equation  [53 

cr  =  y/» 

o-  =  constant 

cr_  ,  cr'  =  Longitudinal  and  circumferential  stress  on 

©  ? 

middle  surface 

"^  =   Azlmuthal  angle 


CJ^      =  Pp 
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NYO-7967 
Buckling  of  Shallow  Spherical  Shells 
under  External  Pressure 

1»   Introduction 

Consider  a  shallow  spherical  shell,  clamped  along 

the  edge  and  under  the  action  of  uniform  external  pressure. 

For  values  of  pressure  increasing  from  zero,  the  shell 

will  at  first  deform  in  a  rather  continuous  manner.  This 

process  persists  \intil  the  external  pressure,  p,  reaches 

a  certain  critical  value,   p^^,   called  the  initial 

cr 

buckling  load.   At  this  value  of  pressure  the  shell  no 
longer  deforms  in  a  continuous  manner  but  jumps  or  snaps 
into  another  non-adjacent  equilibrium  configuration. 
The  pressure  to  which  the  shell  jumps  is  called  the  final 
buckling  load.   The  purpose  of  this  paper  is  to  determine 
the  values  of  p    for  various  shell  configurations. 

A  series  of  experiments  have  been  conducted  by  Kaplan 
and  Pung  (1)  on  clamped,  shallow,  spherical  shells  vinder 
uniform  external  pressure.  Not  only  were  the  initial  and 
final  buckling  loads  determined,  but  measurements  were 
made  of  the  vertical  deflections  at  various  points  of 
the  shell  for  several  different  values  of  the  pressvu?e. 
As  previously  noted  (2,3)>  there  is  an  interesting  inter- 
pretation of  these  experimental  results.   In  Figure  2, 
three  graphs  of  the  vertical  deflection  v  versus  the 
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polar  angle  '5  are  shown.   For  low  values  of  p   (a 
parameter  that  characterizes  the  geometry  of  the  shell, 
see  equation  [53  and  the  text  below  this  equation  for 
the  definition  of  p),  the  shell  deformed  as  in  Figure  2a, 
with  a  maximum  deflection  at  the  center  9=0.   This 
will  be  called  mode  of  deformation  I.   For  p   in  the 
range   20  <  p  <  SSt      the  shell  deformed  in  mode  of  deforma- 
tion II,  as  shown  in  Figure  2b.   This  mode  has  the  maximum 
deflection  point  between  the  center  and  the  edge  of  the 
shell.   For  values  of  p  greater  than  SS,   mode  III  was 
observed,  see  Figure  2c.   In  this  mode,  the  maximum 
deflection  is  again  at  the  center.   Two  important  facts 
about  this  behavior  are:   there  exist  certain  transitional 
values  of  p,   p   equals  approximately  20  and  55,  for 
which  the  mode  of  deformation  changes;  and,  associated 
with  a  change  in  mode,  there  is  a  change  in  sign  of  the 
curvature  at  ©  =  0  for  the  graphs  of  deflection  versus 
5.   That  is,  in  Figure  2,  the  curvatures  at   »  =  0  for 
modes  I  and  III  are  positive,  while  for  mode  II  the 
ctirvature  is  negative.  Therefore  there  should  be  values 
of  p  for  which  the  curvature  at  the  center  vanishes. 
These   p's   are  defined  as  the  transitional  values, 
p^  (i  =  1,  2,  3,  ... ). 

When  the  experimental  initial  buckling  loads  are 
plotted  against  p  rather  than  the  geometric  parameter 
X  of  Kaplan  and  Fung  (1)   (p  =  {2)"'^'^   X^ )  ,   it  is 
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found  (2)  that  there  is  an  oscillating  or  peaking  behavior 
of  these  buckling  loads  as  a  function  of   p,  see  Figure  3. 
In  Figure  3  P    is  the  critical  value  of  the  loading 
parameter,  related  to  the  physical  pressure,   p,   by 
equation  [i+l.   The  circles  are  the  experimental  values, 
while  the  solid  curve  indicates  the  trend  in  the  data. 
It  is  important  to  notice  that  the  peaks  in  the  initial 
buckling  loads  occur  for  those  values  of   p  at  which 
there  is  change  in  mode  of  deformation,  i.e.,  for  the 
transitional  values   p.  • 

Several  approximate  solutions  have  been  given  for 
the  non-linear  boundary  value  problem  governing  the 
deformations  of  a  clamped,  shallow,  spherical  shell.   The 
results  for  a  perturbation  type  of  solution,  with  the  non- 
dimensionalized  radial  displacement  at  some  point  as 
perturbation  parameter  (i|),  and  a  power  series  solution 
(3),  are  also  shown  in  , Figure  3.   The  initial  buckling 
loads  predicted  by  the  perturbation  solution  and  illus- 
trated by  the  dashed  line  are  in  agreement  with  the 
experiments  only  in  the  neighborhood  of  p  =  20.  For 
increasing  values  of  p,   the  agreement  becomes  poorer. 
Furthermore,  the  perturbation  solution  fails  to  reveal 
the  peaking  behavior  of  the  buckling  loads  as  a  function 

of  p. 

The  results  from  the  power  series  solutions  are 
indicated  by  the  short  vertical  lines  in  Figure  3«  The 
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value  of  the  buckling  load  for  a  prescribed   p   lies 
somewhere  in  the  open  interval  spanned  by  the  vertical 
line.   It  can  be  observed  that  the  power  series  solution 
yields  buckling  loads  which  are  in  good  quantitative  agree- 
ment with  experiments.   In  addition  the  peaking  trend  in 
the  buckling  loads  is  also  clearly  evident  from  this  solu- 
tion.  The  disadvantage  of  the  power  series  approach  is 
that  it  yields  discrete  results  rather  than  a  continuous 
variation  with   p.   That  is,  buckling  loads  can  be  computed 
for  only  one  value  of   p   at  a  time.   Furthermore,  for 
larger   p   and  P,   the  niomber  of  terms  in  the  series  and 
the  number  of  iterations  required  for  convergence  of  the 
numerical  procedure  become  quite  large.   The  calculations 
thus  become  rather  cumbersome. 

In  this  paper  a  simplified  approach  to  predicting 
initial  buckling  loads  is  presented.   This  approach  is 
based  upon  the  solutions  to  two  linearized  versions  of 
the  original  non-linear  problem.   One  of  the  linearized 
versions  is  an  eigenvalue  problem  while  the  other  corre- 
sponds to  the  linear  elasticity  problem  for  the  bending 
of  a  shallow  spherical  shell.   Although  the  method  is 
rather  crude,  it  yields,  in  a  simple  manner,  results 
which  are  in  excellent  agreement  with  the  experiments. 
Finally,  a  semi-empirical  formula  is  given  for  determin- 
ing Initial  buckling  loads. 
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2«   The  Nonlinear  Problem 

The  definitions  of  geometrical  quantities  associated 

with  the  spherical  cap  are  shovm  in  Figure  1.   The  polar 

angle  "S,      and  the  angle   a   that  the  tangent  plane  to 

the  deformed  siorface  makes  with  the  horizontal 

^"^  ~  R  ~Z  ■*■©)»   may  be  normalized  as  follows: 
"  d© 

[1]  ©  =  _A.©   ;  a  =    J\.a 

where  -A.  is  the  semi-angle  of  opening  of  the  shell.   If 
a  stress  function  y  is  introduced  such  that 

[2]    0-^ ^—   -  ^   ,   o-^ -—  -  -X 

o        o 

where   (T   and  cr   are  the  longitudinal  and  circumferen- 

tial  membrane  stresses,  then  the  nonlinear  differential 
equations  for  the  rotationally  symmetric  deformations  of 
a  spherical  cap  are  as  given  in  (3) 

[3a, b]    La  =  p(aY  +  P©^ )  }      Ly  =  p(©^  -  a^) 

L   is   the   linear  differential   operator, 

^   '  ®   d©   ©  I©    ^ 
and      P     and     p,      the  loading   and  geometrical   parameters, 
are   defined  by   the   following  relations, 
[h]  P  =  2Ek^/2    (hj2   p 

[5]  p   =   {k)-l/2^2  R 

where 
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k  =  — ^-^ 

3(l-v2) 
and  p  is  the  uniform  external  pressure.  The  geometrical 
parameter  p  is  also  proportional  to  the  ratio  of  the 
shell  rise  and  the  thickness. 

Associated  with  the  differential  equations  [3a, b] 
are  the  following  boundary  conditions:   At  the  center, 
©  =  0,  there  is  rotational  symmetryj  while  at  the  edge, 
9=1,  the  shell  is  clamped  and  the  displacement  is  zero. 
Mathematically  this  may  be  expressed  (see  reference  (3)) 
as 

[6a]  a(0)  =  0   }   y(0)  =  0 

[6b]  a(l)  =1   i        [||  -  VY]e=i  =  0 

Thus,  the  complete  fonnvilation  of  the  nonlinear 
boundary  value  problem  consists  of  the  differential 
equations  [3a,b]  and  the  boundary  conditions  [6a,b]. 
In  reference  (3)  an  approximate  solution  based  upon  a 
power  series  expansion  of   a  and  y   i^i  terms  of  9, 
has  been  presented  for  a  large  range  of  P  and  p  • 
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3»   The  Linear  Problems 

Rather  than  attempt  a  direct  solution  to  the  nonlinear 
boundary  value  problem,  the  initial  buckling  loads   P^^ 
are  determined  from  the  solutions  to  two  linearized  problems 
associated  with  equations  [3a,b]  and  [6a,b].  The  first 
of  these  problems,  which  is  a  linear  eigenvalue  problem, 
will  be  called  Problem  A;  the  second  one,  which  is  the 
linear  elasticity  problem  for  the  bending  of  a  shallow 
spherical  shell,  will  be  called  Problem  B. 

Problem  A: 

To  obtain  the  differential  equations  for  Problem  A, 
it  is  convenient  to  make  the  following  change  of  dependent 
variables: 
[7]  a=y  +  9       Y  =  0"© 

so  that   y  is  proportional  to  the  slope.   Then  the 
equilibrixxm  equation,  [3a],  becomes 

[8]  Ly  =  ©p  (y<r  +  o-©  +  P©)  . 

It  will  now  be   asstimed  that      cr      may  be  written   as 

[9]  o-    =  <r°  +  5<r  (©) 

where  a  °   is  a  constant,  and   5o-   is  small  compared 
to  0-°,   This  is  equivalent  to  assuming,  by  equations 
[2]  and  [?],  that  the  stress  is  constant  correct  to 
higher  order  terms.   From  the  results  of  the  numerical 
solution  of  (3),  this  is  not  an  unreasonable  assumption. 
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By  substituting  [9]  Into  [8]  and  neglecting  higher  order 
terms,  there  results 

d©  © 

The  second  of  the  boundary  conditions  of  [6b]  is 
replaced  by  the  condition  that  the  total  vertical  force 
along  the  edge  9  =  1  is  equal  to  the  vertical  force 
produced  along  the  edge  by  the  external  pressure.   This 
condition  yields  the  result  that  at  ©  =  1, 

[11]  0-°  =  -  P 

Setting 

[12]  cO^  =  Pp 

and  by  substitution  of  [11]  and  [12],  equation  [10] 
becomes 

d©^  © 

The  linear  boundary  value  problem  formed  by  equation 
[13]  and  the  boundary  conditions, 

[lU]  y(o)  =  yd)  =  0 

obtained  from  [6a,b]  with  the  aid  of  [?],  provides  a 
rough  approximation  to  the  original  nonlinear  problem 
[3a, b]  and  [6a, b]. 

The  assumptions  that  have  been  made  in  deriving 
this  approximating  problem  are  based  on  a  consideration 
of  an  analogous  behavior  between  the  spherical  shell 
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under  external  pressure  and  the  beam-column.  The  analogy- 
is  not  complete  In  the  sense  that  in  the  beam-coliunn, 
the  lateral  and  end  loads  may  be  varied  independantly. 
While  in  the  shell,  the  "edge  load"  at   6=1   is 
dependent  on  the  external  pressure   P,   Essentially, 
the  assumption  [9]  and  the  condition  [11]  replace  the 
clamped  spherical  shell  under  external  pressure,  in  a 
first  approximation,  by  a  shell  that  is  clamped  along 
the  edge  ©  =  1,  and  loaded  by  a  uniform  force  propor- 
tional to  (T  •   This  shell  deforms  by  contracting  into 
a  smaller  spherical  segment  (cr  =  const,,  y  =  0)  until 

some  critical  value,   P  _,   of  <r  °   is  reached.   At 

'   or 

this  value  of  cr    the  shell  may  buckle  into  non-spherical 
states. 

This  eigenvalue  problem,  defined  by  equations  [133 
and  [II4.],  has  non-trivial  solutions 

[15]  y  =  C  J^icj  ©)   , 

for  a  discrete  set  of  values,  u   (n  =  1,  2,  ••,),   of 
U)   satisfying  the  equation 

[16]  •^l^^n^  "  ° 

where  J^{u)g)    are  the  Bessel  functions  of  the  first 
kind  of  order  one.   The  trivial  solution  satisfies  the 
compatibility  equation  [3b]  under  assumption  [9]« 
Therefore,  from  the  definition  [12],  the  initial  buckling 
loads  may  be  expressed  as  a  function  of  p  through, 
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[17]  ^cr  =  "T  '  ^  =  1»  2,  3,  ... 

In  the   P    ,   p  plane  equation  [l?]  yields  a  family  of 
equilateral  hyperbolas.   Each  member  of  the  family  repre- 
sents a  different  mode  of  deformation.   The  problem  now 
Is  to  determine  those  values  of  p  for  which  the  solution 
moves  from  one  member  of  the  family  to  another.  I.e.,  for 
which  values  of   p  there  is  a  transition  from  the  ^    _-, 
hyperbola  to  the   ^    hyperbola.   Since  these   p's   repre- 
sent the  point  of  change  from  one  mode  to  another,  they 

-v. 

are  the  transitional  values   p.  •   They  will  be  determined 
from  the  solutions  of  Problem  B. 
Problem  B: 

This  is  the  linear  elasticity  problem  for  the  bend- 
ing of  a  shallow  spherical  shell  under  uniform  external 
pressure.   By  considering  infinitesimal  displacements 
from  equilibrium,  neglecting  higher  order  quantities  in 
equation  [3a,b],  and  then  substituting 

a  =  y  +  9   and   Y  =  P  -  P® 

into  the  resulting  differential  equations  and  the  boundary 
conditions  [6a, b],  the  following  boundary  value  problem 
of  the  linear  bending  theory  is  obtained: 

tlQ^'^^    i©  I  §9  ®7  -  PP  =  0  i  §5  11^  ^P  .  2  py  =  0 

[19a]  y(0)  =  0    i    P(0)  =  0 

[19b] 


yd)  =0   ;    If-  -  vP     =  (1-v)  P 

L      -i©=l 
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The  solution  to  this  problem  y  =  yc(^»P»p)  inay  be 
obtained  in  terms  of  the  Kelvin  functions  ber  ©   and 
bei  9   (ber  £,   =  ReJ^d^/^^),  bei  ?  =  ImJ^d^^^C ) ).   It 
depends  parametrically  on  P   and  p.   The  mode  of 
deformation  is  independent  of   P,   and  varies  only  with 
the  geometrical  parameter  p.   The   pT   have  previously 
been  defined  as  those  values  of   p   for  which  the  curva- 
ture of  the  deformed  state  vanishes  at   9=0.   Since 
y  is  proportional  to  the  slope,  the  p.   are  the  roots 
of 


[20]  dyg 

d9 


=  0   . 
9=0 


For  the  boundary  conditions  [19a, b]  the  condition  [20] 
reduces  to 

[21]  ber'  Z^  =  0 

where 

[22]  Z^  =.2^/^  (pp^/^   . 

The  first  three  roots  of  [21],  (5),  yield  the  transi- 
tional values 

[23]     p'l  =  25.8   ,   p*  =  76.2   ,   P3  =  158.5  . 

The  two  linear  problems  A   and  B   are  independent 
in  the  sense  that  they  do  not  supply  contradictory 
information  in  the  same  area.   Problem  A  was  formulated 
30  as  to  obtain  approximate  values  of  the  buckling  loads 
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considered  as  a  bifxircation  point  in  the  solution.   It 
therefore  could  not  be  expected  to  yield  accurate  Informa- 
tion concerning  deformations.   On  the  other  hand,  Problem 
B  should  yield  reasonable  results  concerning  modes  of 
deformation  and  cannot  be  expected  to  supply  results 
with  respect  to  the  bifurcation  phenpmenon.   Combining 
the  solutions  to  Problems  A  and  E,  equations  [17]  and 
[23],  the  initial  buckling  loads  for  a  clamped  spherical 
cap  may  be  evaluated  from  the  following  formula: 


r     2 
(J 


[2i|] 


where  (5) 


cr 


^/p  for  0  <  p  <  25.8 
|/p  for  25.8  <  p  <  78.2 
CJ^/p  for  78.2  <  p  <  158.5 


U 


[2i^a] 


O^  =  3.832 


,   £v>2  =  7.016 


CJ 


3  =  10.171+ 
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[|.>  Comparison  with  Previous  Results  and  Experiments 

The  values  of  the  Initial  buckling  loads,  equation 
[2lj.],  obtained  from  the  approximating  linear  theories 
are  indicated  by  the  solid  curve  in  Figure  l\..      For 
comparison,  the  results  of  the  perturbation  and  power 
series  solutions  as  well  as  the  experimental  results  of 
Kaplan  and  Pung  are  also  shown  in  Figure  i;.   It  can  be 
observed  that  the  simplified  method  given  in  this  paper 
yields  results  which  are,  especially  in  comparison  with 
previous  solutions,  in  excellent  agreement  with  the 
experiments.   The  major  point  of  disagreement  with  the 
experiments  are  the  numerical  values  of  p,  •  The  values 
computed  by  this  theory  are  greater  than  those  observed 
experimentally*   Two  possible  explanations  of  this  are: 

1)  The   p.  ,   equation  [23],  are  computed  for  a 
shell  that  is  clamped  on  the  boundary,  see  [19a,b], 
In  a  physical  experiment  it  is  impossible  to  provide 
exact  clamping  on  the  edge;  the  actual  condition  on  the 
edge  is  somewhere  between  clamped  and  simply  supported, 

2)  The  differential  equations  [l8a,b]  which  were 

•it 
employed  in  computing  the   p.   are  based  on  the  assump- 
tion of  a  shallow  shell.   The  resulting  p.   [23]  may 
be  in  error  for  the  larger  p,   since  these  values  of 
p  correspond  to  shells  that  are  non-shallow. 

To  illustrate  point  1),  consider  Problem  B  for  a 
simply  supported  edge.   With  the  boundary  conditions 
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[19a,b]  replaced  by  those  of  simple  support,  it  can  be 
shovm  that  the  equation  corresponding  to  [21]  whose  roots 
yield  the  p'.      for  this  case  is: 

[25]   v/2  Z^  bei  Z^  +  (l-v)(ber^  Z^  +  bei^  Z^)  =  0   , 

where   Z.   is  defined  in  equation  [22].   The  first  two 
roots  of  this  equation  are 

[26]         Pi  =  17    ,   p|  =  62.3   . 


Thus  by  considering  a  non-clamped  edge,  the  values 
of  p.   are  decreased.  It  may  be  noticed  from  Figure  Ij. 
that  the  experimental   p,   is  bracketed  between  the 
clamped  and  simply  supported  p,  •   This  justifies  the 
conjecture  concerning  the  effect  of  variation  of  boundary 
conditions  between  theory  and  experiment.  The  further 
disagreement  between  theory  and  experiment  in  pp   is 
probably  caused  by  the  non- shallowness  of  the  shell  for 
this  value  of  p. 

On  the  basis  of  this  evidence,  the  following  semi- 
empirical  formula  is  offered  in  place  of  [2l|.]  for  comput- 
ing the  initial  buckling  load  for  clamped,  shallow 
spherical  shells: 

CJ  J/p  for   0  <  p  <  20 

[27]       Pjjj,  =  -(  CJ^/p  for  20  <  p  <  56 

(J^/p  for  56  <  p  <  N 
where,  U)        are  given  by  equation  [2]4.a],  and  N   is  to 
be  determined  from  experiments. 

-  18  - 


BIBLIOGRAPHY 


!•      Kaplan,   A.,    and  Fving,    Y.    C,    A  Non-linear  Theory 
of  Bendinp;   and  Buckling  of  Thin  Elastic   Shallow 
Spherical   Shells.    NACA  TN   3212,    193U» 

2.  Relss,  E,  L.,  On  the  Non-linear  Buckling  of  Shallow 
Spherical  Domes,  Jour.  Aero»  Sol,,  Vol.  23^  No«  10, 
pp.  973-975,  October,  1956. 

3»   Relss,  E.  L.J  Greenberg,  H.  J,,  and  Keller,  H.  B., 
Non-linear  Deflections  of  Shallow  Spherical  Shells, 
Report  NYO-7695  of  AEG  Computing  Facility,  Insti- 
tute  of  Mathematical  Sciences,  New  York  University, 
published  December  1956.   Paper  presented  before 
25  Annual  Meeting  of  Inst,  of  Aero.  Sci.,  Reprint 
No.  663,  January  1957. 

k..      Archer,  R,  R.,  Stability  Limits  for  a  Clamped 

Spherical  Shell  under  Uniform  Pressvire,  Doctoral 
Thesis,  Mass.  Inst,  of  Tech.,  May,  1956. 

5«   McLachlan,  N.  W.,  Bessel  Functions  for  Engineers, 
Oxford  Univ.  Press,  London,  193i4-» 


-  19  - 


h=  ^thickness 


SHELL  GEOMETRY 


FIG.  I 


20<o<55 


V       \  55<yO<? 


FiG.    2 

MODES  OF  DEFORMATION 


COMPARISON  OF  PREVIOUS   THEORETICAL 
AND     EXPERIMENTAL    INITIAL   BUCKLING 

LOADS 


CLASSICAL  THEORY  FOR    SPHERE 


PCR 


oEXPI     REF  I 

—  TREND  OF    EXPTL 
DATA 

PERT  SOL. 

REF    4 

POWER  SERIES  SOL. 
REF  3 


100 


P 


FIG   3 


COMPARISON  OF   INITIAL 
BUCKLING  LOADS 


CLASSICAL  THEORY  FOR  SPHERE 


p 
XR 


FIG.  4 


NEW  YORK  uNtveRSin 

INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

25  WaverV  Place,  New  York   3,  N.  Y, 


NYU 

NYO- 

7967       Reiss. 

Bucklino;  of  sh«nr.T., 


NYU 

NYU 
NYO- 

AlIthcJr 


^eiss, 


-^^ic^llng^of  shallow 
^4fpherical   sE^11^~tIHI^ 


DATE   DUE 


BORROWERS    NAME 


ROOM 

NUMBER 


N.  Y.  U.  Institute  of 
Mathematical  Sciences 

25  Waverly  Place 
New  York  3,  N.  Y. 


